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The tilt-over mode in a precessing triaxial ellipsoid is studied theoretically and numerically. In- 
viscid and viscous analytical models previously developed for the spheroidal geometry by Poincare 
[Bull. Astr. 27, 321 (1910)] and Busse [J. Fluid Mech., 33, 739 (1968)] are extended to this 
more complex geometry, which corresponds to a tidally deformed spinning astrophysical body. As 
confirmed by three-dimensional numerical simulations, the proposed analytical model provides an 
accurate description of the stationary flow in an arbitrary triaxial ellipsoid, until the appearance at 
more vigorous forcing of time dependent flows driven by tidal and /or precessional instabilities. 
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I. INTRODUCTION 

The flow of a rotating viscous incompressible homogeneous fluid in a precessing container has been studied for over 
one century because of its multiple applications, such as the motions in planetary liquid cores and the generation of 
planetary magnetic fields [e.g. UQ- In the spheroidal geometry, the early work of Q demonstrated that the flow of 
an inviscid fluid has a uniform vorticity. Later, viscous effects have been taken into account as a correction to the 
inviscid modes, in considering carefully the critical regions of the Ekman layer Indeed, the Poincare solution is 

modified by the apparition of boundary layers, and some strong internal shear layers are also created in the bulk of 
the flow, which do not disappear in the limit of vanishing viscosity @. Besides, at high enough precession rates, these 
shear layers may become unstable 0-0]; an d m a second transition, the entire flow becomes turbulent: this is the 
. precession instability. These experimental works have been completed by numerical studies in cylindrical, spherical 
Oh' and spheroidal geometries [IB; 14 L in particular for studying kinematic dynamo models in a spheroidal galaxy [l5| or 
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for geophysical applications 16-20}. 

However, in natural systems, both the rotation and the gravitational tides deform the celestial body into a triaxial 
\ ellipsoid, where the so-called elliptical (or tidal) instability may take place (see [2l| - [25| for details on this instability and 
its geo- and astro-physical applications). The elliptical instability takes place in any rotating fluid whose streamlines 
are elliptically deformed (see e.g. [26[ or (2?J )■ It comes from a parametric resonance of two inertial waves of the 
rotating fluid with the tidal (or elliptical) deformation of azimutal wavenumber m = 2 [25} . Similarly, it has been 
suggested that the precession instability comes from the parametric resonance of two inertial waves with the forcing 
related to the precession of azimutal wavenumber m = 1, which comes from the deviations of the laminar tilt-over 
base flow from a pure solid body rotation [28l - [3fj| . But it has also been suggested that the precession instability is 
related to a shear instability of the zonal flows that appear in a precessing container [e.g. [1^]. Clearly, the precise 
. . ■ origin of the precession instability is still under debate, and is beyond the point of the present work. But since tides 
and precession are simultaneously present in natural systems, it seems necessary to study their reciprocal influence, 
in presence or not of instabilities. The full problem is rather complex and involves three different rotating frames: the 
precessing frame, with a period T p w 26000 years for the Earth, the frame of the tidal bulge, with a period around 
Td « 27 days for the Earth, and the container or 'mantle' frame, with a period T s 23.93 hours for the Earth. As 
a first step towards the full study of the interaction between the elliptical instability and the precession, we consider 
the particular case where the triaxial ellipsoid is fixed in the precessing frame (T<j = T p ), which allows the theoretical 
approach to be analytically tractable. We show in figure [T] a sketch of this configuration. 

The paper is organized as follow. In section UH Poincare and Busse analytical models are extended to precessing 
triaxial ellipsoids. Then in section Hill our analysis is validated by comparison with a numerical simulation. 
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FIG. 1: Sketch of the problem under consideration. A hollow solid but deformable spheroid (i.e. 'the mantle') is filled with liquid 
and set in rotation at a constant angular velocity Q versus its axis (Ox-j). The spheroid axis is tilted at the precession angle 9 
and fixed on a rotating table, which rotates at the precession rate £l p . Two fixed rollers aligned with (Oxs) then transform the 
spheroid into a triaxial ellipsoid by compression along the axis (Oxz), perpendicular to the rotation axis. 



II. ANALYTICAL SOLUTION OF THE FLOW IN A PROCESSING TRIAXIAL ELLIPSOID 

In this section, we consider firstly the case of an inviscid fluid and extend the Poincare model to a processing triaxial 
ellipsoid, which allows us to obtain explicit analytical solutions. We then tackle the viscous case in extending the 
Busse model, following the method of As sketched in figure [TJ we consider the rotating flow inside a precessing 
triaxial ellipsoidal container of principal axes (01,02,03). We define (Oxi, Ox2, Ox%) in the frame of the tidal bulge 
which is also the precessing frame, such that Oxi is along the principal axis Oj of the ellipsoidal container and (Ox^) 
is the mantle rotation axis. We note f2 the imposed mantle angular velocity and use f2 _1 as a timescale. We also 
introduce the mean equatorial radius R eq — (ai + a2)/2, which is used as a lengthscale. Consequently, the problem 
is fully described by six dimensionless numbers: the eccentricity e = (a\ — a\)/(a\ + a|), the aspect ratio 03/01, 
the Ekman number E = v/Q,B^ q , where v is the kinematic viscosity of the fluid, and the three components of the 
dimensionless precession vector ti p in the inertial frame of reference, i.e. the angle 9 between H p and e X3 , the angle 
82 between ft p and e Xl and the angular precession rate Cl p , which is positive for prograde precession and negative for 
retrograde precession. 



A. Inviscid Poincare tilt-over mode for a triaxial ellipsoid 



We look for a base flow solution of the Euler equations, i.e. 

du 

— + u • Vu = —Vp - 2J1 P x it, 



(1) 



Vu = 0, 



(2) 



inside the triaxial ellipsoid in the precessing frame. In this frame, using the ellipsoid equation 




fce(i,2,3) 



and we write the 



(3) 




(4) 
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As suggested in Q, we focus on the so-called 'simple motions' such that the velocity U(Ui, U%, U3) is described as 
linear combinations of the coordinates. This hypothesis leads to a solid body rotation in the sphere, i.e. U' = ut' X r', 
where u/(u4,u; 2 , oj' 3 ) depends on time a priori. Equation |4] then gives the following velocity field in the ellipsoid: 



u - (^fc)fce(i,2,3) = ( — — ^fc+i^fc-i — ^fc-i^fc+i J 

v ; \ak-i afc+i / 



(5) 

fee(i,2,3) 



where permutations k G (1, 2, 3) are used. 

Now, we have to find the 'simple motions' solution of the Euler equations for the rotational part of the flow, taking 
into account the no-penetration boundary conditions for its irrotational part, which leads to the so-called Poincare 



flow. We follow the method of [3l| rather than the lagrangian method of Poincare which is more laborious. We deduce 



the rotation rate vector ui from the velocity field (JS|): 

- = \ V x U = K) fce(1 , 2 , 3) = U( ^ + ^) <J h ) - (6) 

2 2 \\a k+1 a k -ij / k e(l,2,3) 

Note that u>' is independent of the space coordinates so that u> is uniform. Taking the rotational of ([T]), the inviscid 
equation for u) gives the following three scalar equations: 

—j^ + (a/c+i,fc — aik-i,k) LOk+i = &k-i,k &p,k-i Wfc+i — oik+i,k Qp,k+i (7) 

2 2 

for permutations k € (1, 2, 3), with the coefficients a.ij = — 2 +2 = 2 ^ J '^ 2 and the different ellipticities r/ij = — J? 3 of 

the container. The equations for a spheroidal geometry arc recovered with a\ = a 2 . 

We focus on the stationary solutions of the problem. In this particular case, we solve the system ([7]) analytically, 
which gives: 

w 3 aj + al 

Wi — S2 Pl i, (8) 

Pi 2 7l 

u 2 — -3 " P ,2, (.yj 

P12 72 

where = 2 ' a — 3 — /3ji and 7, = jf^{aj — oj) + 2 f2 Pi 3 ai a 2 . The corresponding velocity field writes: 

E7i = 01 — + , (10) 

P12 \ a 2 72 / 

rr w 3 (xi 2 a 2 ftp.l^N 

^ 2 = a 2 — 2 , (11) 

P12 \ a i 7i / 

t/3 = «3 -5— ■ (12) 

P12 V 71 72 / 

Note that the choice of uj 3 is arbitrary here. Actually, lo 3 is determined by the boundary (Ekman) layer and thus 
has to be determined by a viscous study, following for instance the method of ||, as described in the next section. 
Note also that this velocity field is divergent for 71 = or 7 2 = 0: the inviscid study gives two resonances, which 
correspond to resonance between the frequencies of respectively the precessional forcing and the tilt-over (see i for 
details). These linear resonances are reached for two specific precession rates (depending on the aspect ratio) which 
are 



a 



^.3 = ^3 (13) 



for i £ (1,2). It is clear that oblate ellipsoids (ai,a 2 > 03) have their resonance in the retrograde regime (i.e. in 
the range f2 Pi 3 < 0) whereas prolate ellipsoids (ai,a 2 > 03) have their resonance in the prograde regime. Finally, 
compared to the spheroidal case, an important result here is the apparition of the second resonance created by the 
equatorial ellipticity of the container. 
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B. Viscous study of the tilt-over mode in a triaxial ellipsoid 



Following |5|, it is possible to take into account the viscosity in the study of the flow in a precessing triaxial ellipsoid. 
Here, we focus on the equivalent method of Q, based on the equilibrium between the inertial torque Tj, the pressure 
torque T^, and the viscous torque T v . Keeping the leading terms for these three torques, the general torque balance 
given in [5| for a steady rotating flow q = u> x r in the precessing frame within a volume V with a surface £ writes: 



2 / r x (fl p x q)dV = - f pr x nd£ + E rx V 2 qdV, 
Jv /s Jv 



where r is the position vector and n is the unit vector normal to £ pointing outward [see[ 
these terms for a triaxial ellipsoid. 

At first order, the pressure gradient equilibrates the centrifugal force, which gives: 



(14) 



Now, we have to calculate 



+ wf_ t ) x\ - Wi Uj X, Xj 



(15) 



Hence in the limit of small ellipticities, i.e. at first order in 77 = 1 — 03/01 and 772 = 1 — 02/01, the pressure torque 
writes: 



(77 - i] 2 ) lo 2 w 3 

T p = - (f> prxn dS = (r p ,fc) fce ( li2 ,3) = ^ I -77 wi w 3 

772 Wi w 2 



(16) 



where I is the moment of inertia in the spherical approximation. Note that the expression of 
spheroid (i.e. ai = 02 hence 772 =0). 

After little algebra, the precessional torque simply writes: 



= 7fi p x u. 



is recovered for the 



(17) 



Finally, the equations (TT6l) and (fl~7| gi 



u> = and 1^ • u) — 0, and thus with the equation (fl4| F v ■ u> = 



Consequently, the fluid being in a stationary state, there is no differential rotation along u>: 



w • e 3 



(18) 



Indeed, in the rotating frame of the fluid, the angular rate of the container along oj is null: the only differential 
rotation between the fluid and the container is in the equatorial plane such that no spin-up process occurs. Note 
that this can also be recovered with a boundary layer analysis, in the same way as which shows that the volumic 
Ekman pumping is solution of the inviscid bulk equations provided that the so-called solvability condition (TT51) is 
verified. According to this equation, also named the no spin- up condition in Q , the only relative motion between the 
interior and the boundary is a rotation given by u) eq = u> — ft. Then, we need to calculate the viscous torque due 
this equatorial differential rotation. This calculation relies on the fact that in order to maintain the basic stationary 
state, the torque supplied to the fluid to counter balance the viscous friction is given by the decay rate that would 
occur if at a given time the precession is turned off. This decay rate is simply given by the rate at which energy is 
dissipated by the Poincare mode in a free system at t = 0. This rate is given by the Greenspan's theory Q, valid in 
the frame rotating with the fluid. Thus, this linear solution for the viscous decay of the spin-over mode in a rotating 
fluid leads to introduce a new Ekman number Ef — E/u and a new unit of time t = iw, scaled with the fluid rotation 
rate ui. According to @, the time evolution of uj eq in the non-rotating frame is: 



cos (\i i \f~Epj w e q(0) — sin ^Ai t \[E~f 



iv x iv eq (0) 



(19) 



with A r = —2.62 and Ai = 0.259. Note that strictly speaking in triaxial ellipsoids, the ellipticity modify the growth 
rate and eigenfrequency of inertial modes. However, in the limit of small ellipticities we consider here, this modification 
can be neglected. 
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FIG. 2: Theoretical amplitude of tilt-over mode angular rotation rate u) for a spheroid (e = 0) and for a slightly deformed triaxial 
ellipsoid (e = 0.03). The other parameters used in this figure are those of J^j: the rotation rate of the container fl — 207 rpm 
and a\ = 0.125 m (such that the Ekman number is E = 3 ■ 10~ 5 ), 8 — 9°, a^/ai = 0.96 and 62 = 0°. The spheroidal case 
corresponds to the plot v — 10~ 5 m 2 / s in figure 3 of f^j. Note that even a very small tidal deformation significantly changes 
the obtained solution. 



Using this equation, reintroducing the variables E and t, the equatorial viscous torque is: 



/ du) eq 
V dt 



IVloE I A r u>2 — K uji/u I (20) 
t=o \ \ r (u 3 - 1) 



Then, the torque balance given by f| 14[) projected onto the rotation axis of the fluid oj (the no spin-up condition 
T8| ). as well as onto the principal axes e Xl and e X3 yields the following system of equations: 

wj+w2=w 3 (I-W3) (21) 



flp.2 uj 3 - fi Pi3 u 2 = {rj - 772) U>2 w 3 + 



\ 1/4 
A r UJ\ U)^ 



Xi 



UJ 2 

"774 



(22) 



u 2 — flp.2 ui = ri 2 wi 0J2 — A,, cul^ 4 (1 — u> 3 ) \fE. (23) 

The supplementary terms compared to Q or Q do not allow to simplify the system of equations into only one, 
and the full system has to be solved numerically to obtain the rotation axis components of the fluid. This non-linear 
system can be solved in an efficient way with a continuation method (successive perturbations on the a 2 axis) starting 
from the Busse's solution in a spheroid. An example is shown in figurc[2]where the solution in the spheroidal geometry 
(the case v = 10 ~ 5 m 2 /s in figure 3 of @, which gives a Ekman number ofi? = 3-10~ 5 )is compared to a slightly 
deformed triaxial ellipsoid (e = 0.03) with the same ratio a 3 /a\. It can be noticed that even a very small tidal 
deformation e radically changes the obtained solution. 



III. NUMERICAL AND EXPERIMENTAL VALIDATION 



Our purpose here is to validate and test the range of validity of our analytical solution by comparison with numerical 
simulations of the full non-linear Navier-Stokes equations in a precessing ellipsoid. 



6 



A. Numerical resolution 



We consider the rotating flow inside a triaxial ellipsoidal container of principal axes (a\, 0,2, 03), as sketched in figure 
[TJ We work in the precessing frame of reference. Starting from rest, a constant tangential velocity U {x^/a^) 2 
is imposed from time t = all along the outer boundary in each plane of coordinate X3 perpendicular to the rotation 
axis (OX3), where U is the imposed boundary velocity at the equator. We introduce the timescale f2 _1 by writing the 
tangential velocity along the deformed outer boundary at the equator U — QR eq - We then solve the Navier-Stokes 
equations with no-slip boundary conditions, taking into account a Coriolis force associated with the precession O p , 
i.e. in the frame (Oxi, Ox2, Ox^) of the tidal bulge which is also the precessing frame, we solve: 

+ u ■ Vu = -Vp + E A u - 2fl p x u, (24) 



V • u = 0. (25) 

In this work, the range of parameters studied is E > 10~ 3 and e < 0.32. Once a stationary or periodic state is 
reached, we determine the rotation rate in the bulk of the fluid, i.e. outside the viscous boundary layer. To do so, 
we introduce an interior homothetic ellipsoid, in a ratio k, and we define the bulk rotation rate us as the mean value 
of rotation rate over this homothetic ellipsoid. Following [32| , we consider a dimensionless viscous layer thickness of 
8 V fa 5\AB and thus choose for the homothetic ratio k « 1 — i„ « 1 - h\[E. Note that the spheroidal case can be 
efficiently solved by spectral methods (see e.g. [lj| or (20l|). But for the triaxial ellipsoids we are interested in, there is 
no simple symmetry. Our computations are thus performed with a finite element method, which allows us to correctly 
reproduce the geometry and to simply impose the boundary conditions. The solver and the numerical method are 
described in details in I24II. 



B. Experimental set-up 

The experimental set-up has been described in [30l . [33[ for a precessing cylinder and readers should refer to these 
papers for more details. For this paper, the experiment has been slightly modified in order to study the precession of 
a spheroid. This allows to validate the flow at small Ekman numbers, which is impossible numerically. Unfortunately, 
this set-up is limited to a spheroidal geometry {a\ — a?) and is not able to validate the theory for a tri-axial ellipsoid. 
Further modifications to the set-up, such as two rollers compressing the spheroid, would be needed to make the 
equatorial plane elliptical. 

The home-made spheroid has been obtained by assembling two half-spheroidal cavities drilled in solid Plexiglas 
cylinders. The accuracy of the machines (10 microns) ensured that the step between the two parts would be smaller 
than the Ekman layer (of the order of 300 microns at E = 10~ 5 ). This spheroid of equatorial diameter 17 cm and 
aspect ratio 03/01 = 0.85 is filled with water and mounted on a motor which is itself located on a rotating platform. 
The angular velocities of the spheroid and the platform are stable within 0.1% and the precessing angle 8 between 
the two axes was varied from 5° to 15° with an accuracy of 0.1°. 

Pulsed Yag lasers are used to create a luminous sheet perpendicular to the axis of the rotating platform. A 
PIV (Particle Image Velocimetry) camera is located on the rotating platform, aligned with the axis of rotation of 
the spheroid. This allowed to obtain PIV measurements in a plane almost parallel to the equatorial plane (xi,X2) 
and located at a distance x$ = 5 cm above it. Since the Yag lasers are not located on the rotating platform, the 
measurement plane is tilted with an angle with respect to the equatorial plane, which introduces an error of up to 
2cm (for 8 = 15°) in the axial location of the velocity vectors. However, the measured velocity components exactly 
correspond to the x\, x-x components of the velocity because the camera is aligned with the axis of the spheroid. There 
is no distortion of the images at the air-Plexiglas interface (because it is a plane) but there are distortions of the 
images due to the Plexiglas-water spheroidal interface. These deformations are small because the refractive index of 
the Plexiglas and the water are close. They were calculated analytically and checked experimentally using a grid. 
They are located mostly at the boundary of the spheroid and they introduce a maximum error of 15% on the radial 
displacement of the particles. This does not bias the measurements because only the central region of the velocity field 
was used for the treatment of the data. The PIV images were rotated numerically in order to remove the background 
rotation of the flow before being treated by a home-made cross-correlation PIV algorithm. 

In the frame of reference of the spheroid, the 2D velocity field was found to be a nearly uniform translation flow 
whose direction and amplitude vary with the precession frequency Q p . This mean flow is linked to the equatorial 
component of the angular rotation rate (wi,^) which creates a uniform translation flow («3i X3 L02,— 1232 £3 ^1) in 
the plane x% = 5 cm. This flow corresponds to the last terms found in equations (|10p and in the case a% =02. 
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The first terms of these equations correspond to the axial angular rotation (0^3), which in the frame of reference of 
the spheroid is equal to w 3 - 1 in dimensionless form and is thus small outside of the resonances. The measurements 
of the mean velocity and of the mean vorticity thus give in a simple manner the three components of the rotation 
vector u>. Such measurements have been done for the first time in a processing spheroid, and will be compared to 
theoretical and numerical results in the following. 

C. Validation in the spheroidal case 

Before dealing with the problem of precessing triaxial ellipsoids, a first step of this work has been to simultaneously 
check the validity of our numerical tool and the validity of the theoretical development of [B[ over an extended range 
in Ekman number. To do so, we combine numerical and experimental approaches. 

In figure [21 the theory proposed in Q is validated over more than three decades of Ekman numbers for two different 
angles of precession. The interest of the presented results is twofold. Firstly, they validate our numerical model, which 
will now be used to study triaxial ellipsoids. Secondly, this completes the previous validations of Busse's theory [e.g. 
S [HI; in particular for rather large Ekman numbers. Note that, as already noted in [ljj, the theoretical analysis 
requires that the angle between the rotation vector of the container and the rotation vector of the fluid is small 
since otherwise the Ekman boundary layer analysis is no longer valid. This could explain the differences between 
the experimental determinations and the theoretical predictions at the highest precessional forcing angle (9 = 15°). 
A supplementary confirmation of both the numerical model and Busse's theory at such Ekman numbers is given in 
figure [4] (a) where the components of the tilt-over rotation axis are in an excellent agreement for a large range of 
precession rate. At large precession rates, the inviscid results of 0] are recovered. 

D. Numerical simulations of precessing triaxial ellipsoids 

A first series of simulations have been performed for a^/ai = 0.86, 9 = 10°, 92 — 45°, E — 1/600, e = 0.1 and 
various precession rates. Results are presented in figure 2] (b). An excellent agreement is found between the numerical 
simulations and the analytical viscous solution all along the explored range, which demonstrates the validity of our 
extended viscous theory. Note also that the triaxial Poincare inviscid flow is recovered far from the resonances. As 
already noted above, an important feature of the triaxial geometry is the apparition of a second resonance. As already 
noted in the section Hi A[ according to the equation (|13p . the two resonances are in the retrograde regime here because 
ai,d2 > 03- Note also that as already described in the literature (see e.g. Q), the viscosity naturally smoothes the 
inviscid resonance peaks but also modifies their position. 

The tilt-over mode is thus well described by our analytical model at relatively large Ekman number and relatively 
low precession rate and tidal eccentricity. However, two types of instabilities can be expected for more vigorous flows: 
the tidal or elliptical instability at relatively small Ekman and/or large eccentricity e, and the precession instability 
at relatively small Ekman and/or large precession rate. Focusing on the tidal instability, we keep the Ekman number 
above the precession instability threshold, and we increase the eccentricity above the elliptical instability threshold. 

The appearance and form of the tidal instability in the presence of a rotating tidal deformation aligned with the 
rotation axis have been studied elsewhere [23|, Hij . This case corresponds in our notations to the specific configuration 
9 = 0. The most well-known mode is the so called spin-over mode, which appears for low values of tidal bulge rotation. 
As shown in appendix [Vj a theoretical stability analysis demonstrates that this mode may grow upon the Poincare 
flow. Note that the spin-over mode, which is an unstable mode of the tidal instability, and the tilt-over mode, which 
is the base flow of precession, are actually identical, hence should both be described by our stationary analytical 
solution. This is indeed the case, as shown in figure [5] (a) around S7 p = where numerical and analytical solutions are 
in good agreement. Moreover, by comparison with the known results of the elliptical instability in the special case 
9 = (see [23| and for details), we expect the elliptical instability to lead to an unstationary mode for Sl p below 
the resonance band of the spin-over, and to no elliptical instability for Q p above the resonance band of the spin-over. 
These results are also recovered in figure [5] (a) , where oscillatory motions are indeed observed for Ci p 6 [—0.5; —0.2], 
whereas good agreement is found for tt p > between analytical and numerical results. Finally, a systematic study of 
the excited mode frequency as a function of the Ekman number for Q p = —0.4 is shown in figure [5] (b) . As expected 
from the stability analysis presented in appendix [Vj the transition from a stationary flow to an unstationary one is 
found below E w 1/300. 
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FIG. 3: Comparison of the theoretical rotation rate components of Busse's solution for a spheroid (continuous lines) with 
experimental PIV measurements and numerical simulations performed over a large range of Ekman numbers for Sl p = —0.14 
and 03/01 = 0.85. (a) For a precession angle of 5° . (b) For a precession angle of 15° . 
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FIG. 4: TTie rotation rate components of the flow in a spheroid, shown in (a), are compared in figure (b) with those in the 
ellipsoid with the same aspect ratio a^/ai — 0.86 (8 = 10°, 62 = 45° and E — 1/600,). The theoretical rotation rate components 
of the inviscid Poincare solution and of its extension to triaxial ellipsoids are given by the black dashed and dotted lines, 

a 2 1 a 2 

assuming that 103 = /3i2 = ° 2 (see avvendix IT1). The theoretical rotation rate components of the Busse viscous solution 
and of its extension to triaxial ellipsoids are represented by the continuous lines. The two inviscid resonances 71 and 72 are 
represented by the vertical black lines, (a) In the spheroid, the validity of our numerical model is confirmed over a large range 
of precession rate by its good agreement with Busse solution, (b) In the ellipsoid with the same aspect ratio but with e — 0.1, 
our analytical solution is in good agreement with the numerical results. 



IV. CONCLUSION 



This paper presents the first analytical solution and the first numerical simulations of the precessing flow inside 
a triaxial ellipsoid. The extension of Busse viscous solution to this more complex geometry is shown to provide 
an accurate description of the stationary flow, which may be destabilized by oscillating modes of the tidal and/or 
precession instabilities for more vigorous forcing. Clearly, the complete study of instabilities in such a system deserves 
more work and will be the subject of future studies with increased numerical power and a new experimental set-up. 
But the study of the stationary flow performed here already highlights the fact that even a very small tidal deformation 
significantly modifies the resulting flow. This conclusion is directly relevant to the dynamics of planetary cores and 
atmospheres, where precession and tidal deformation simultaneously take place. 
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FIG. 5: We consider the case 8 — 10° , az/ai = 0.86, #2=0 and e = 0.317 (corresponding to (oi +tt2)/2 = az). (a) Evolution 
of the rotation rate components of the flow with the amplitude of precession for a fixed value of the Ekman number E — 1/600. 
Legend is the same as in figure^ (b). Filled symbols correspond to cases where the flow is oscillatory whereas empty symbols 
corresponds to a stationary flow. Note that the sign reversal at Q p = is due to the two possible symmetric orientations of the 
spinover mode of the tidal instability, (b) We now focus on the fi p = —0.4 case where the flow is oscillatory at E — 1/600 as 
shown in (a), and we study the evolution of the pulsation f of the flow with the Ekman number. The flow is stationary for an 
Ekman number larger than E — 1/300, and then oscillates at a fixed pulsation f w 0.67. 



V. APPENDIX: INSTABILITY OF THE INVISCID POINCARE FLOW IN A TRIAXIAL ELLIPSOID 

The work [34| have recently studied the stability of a rotating flow with elliptical streamlines in the particular 
case where the rotation axis executes a constant precessional motion about a perpendicular axis. Using the WKB 
method, (34| achieve to quantify both the influence of this particular precession on the elliptical instability, and the 
influence of the strain on the Coriolis instability for such an unbounded cylindrical base flow. In this appendix, we 
consider the stability of the inviscid Poincare flow in triaxial ellipsoids against small inviscid perturbative rotations. 
The theoretical expression of the growth rate for perturbations that are linear in space variables (i.e. corresponding 
to the classical spin-over mode) can be readily obtained with the same method as in (35[, but with the base flow U 
obtained in section Til Al using an arbitrary value of 0J3, As remind in [25j . the most general perturbation- flow linear 
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in the spatial coordinates is 



Uk -U k = K l+ i(t) —1— Xi-l - Ki-\{t) — — .r, i 



di+l 



(26) 



where the scalar amplitude Ki of these small inviscid perturbative rotations can be written Ki = e ki e at . The total 
flow u have to satisfy the inviscid vorticity equation. Then, U being solution flow at order 0, the solution flow at first 
order in e is solution of 



M k = 



where k = 2 3) and M is a 3 x 3 matrix given by 



(27) 



M 



A12 -A 2 i --B12 a 
Ci B 13 a D12 
-B23 <J C2 D21 



where 



(28) 



At 



2n - ^ 



B, 



di a.j 
H 



a 3 a. 



The solution flow at first order is non-trivial if det(Af) = 0, and the growth rate a is obtained in solving this 
equation. Actually, the 3 rd degree equation det(M) = writes in the simple form a 3 +p a = 0, and the growth rate 
is given by: 



A12B13D21 + A21B23D12 — C1C2B12 



(29) 



-B13-B23-B12 

Note that the dimensionless base flow of [35[ is recovered for the particular value CJ3 = /3i2 . It can also be noticed 
that the extension of the result of [35( to the particular case of a rotating tidal bulge at O p = (0, 0, ttbr) given in [24|, 
is recovered: 



C1C2 
B13B 



(a\ — a§ + 2aia2^tbr){a\ — a\ + 2a\a2^tbr) 



2:; 



{a\ + a|)(a| + a|) 



(30) 



Note finally that this method gives the theoretical inviscid growth rate, which has then to be corrected by a viscous 
surfacic damping term [36j j . where \ is a constant \ ~ 2.62. The threshold is then given explicitly by 



E r = 



1 A12B13D2I + A21B23D12 — C1C2B12 



B13B23B12 



(31) 



For the case presented in figureEl i.e. 9 — 10°, 03 = (ai + a%)/2, 9 2 — and e = 0.317, we predict a critical Ekman 
number for instability around E c « 1/300, and for E = 1/600, we predict an unstable spin-over mode in the range 
ftp G [-0.13; 0.12]. 
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